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ABSTRACT 

We have measured central line strengths for a magnitude-limited sample of early- 
' type galaxies in the Fornax cluster, comprising 11 elliptical (E) and 11 lenticular (SO) 

\ galaxies, more luminous than Mb — —17. When compared with single-burst stellar 

population models we find that the centres of Fornax ellipticals follow a locus of 
' fixed age and have metallicities varying roughly from half solar to two times solar. 

The centres of (lower luminosity) lenticular galaxies, however, exhibit a substantial 
spread to younger luminosity- weighted ages indicating a more extended star formation 
history. 

Galaxies with old stellar populations show tight scaling relations between metal- 
line indices and the central velocity dispersion. Remarkably also the Fe-lines are well 
\ correlated with ctq ■ Our detailed analysis of the stellar populations suggests that these 

scaling relations are driven mostly by metallicity. Galaxies with a young stellar com- 
' ponent do generally deviate from the main relation. In particular the lower luminosity 

, SOs show a large spread. 

■ Our conclusions are based on several age/metallicity diagnostic diagrams in the 

Lick/IDS system comprising established indices such as Mg2 and H/3 as well as new and 

Q^. more sensitive indices such as H7A and Fe3, a combination of three prominent Fe-lines. 

I ' The inferred difference in the age distribution between lenticular and elliptical galaxies 

O I is a robust conclusion as the models generate consistent relative ages using different 

age and metallicity indicators, even though the absolute ages remain uncertain. The 
^ ' absolute age uncertainty is mainly caused by the effects of non-solar abundance ratios 

which are not yet accounted for by the stellar population models. 

Furthermore we find that elliptical galaxies and the bulge of one bright SO are 
overabundant in magnesium, where the most luminous galaxies show the strongest 
, overabundances. The stellar populations of young and faint SOs are consistent with 

■ solar abundance ratios or a weak Mg underabundance. Two of the faintest lenticular 

galaxies in our sample have blue continua and extremely strong Balmer-line absorption 
suggesting star formation < 2 Gyrs ago. 

Key vifords: galaxies: abundances - galaxies: clusters: individual: Fornax - galaxies: 
formation - galaxies: elliptical and lenticular - galaxies: kinematics and dynamics 



1 INTRODUCTION 



Great efforts have been made in the last few years to develop 



evolutionary stellar population synthesis models ( Bruzual & 


Chariot 1993; Worthey 1994 


; Weiss, Peletier & Matteucci 


1995; Vazdekis et al. 1996; 


Kodama & Arimoto 1997) in 



order to analyze the integrated light of galaxies and de- 
rive estimates of their mean ages and metal abundances. 
One of the main obstacles in the interpretation has been 
the age/metallicity degeneracy in old stellar populations. 
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As pointed out by Worthey (1994) the integrated spectral 
energy distribution (SED) of an old (> 2 Gyrs) stellar pop- 
ulation looks almost identical when the age is doubled and 
total metallicity reduced by a factor of three at the same 
time. Therefore two galaxies with almost identical broad- 
band colours can have significantly different ages and metal- 
licities. In the optical wavelength range, only a few narrow 
band absorption line-strength indices and the 4000 A break 
(see also Gorgas et al. 1999) have so far been identified which 
can break this degeneracy. 

One of the most successful and widely used methods 
for measuring the strength of age/metallicity discriminating 
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abs orption features is the L i ck/IDS system ( Burstein et al 
198 4 hf^orthey et al. l"994| ; |Traser et al . 1998|) which has 



been u sed by m a ny authors (^ onzalez^^993 ; Davies, Sadler 
Pelel |ier 199j; [Fisher, Frarix fc Illingworth 1995|; [Fisher 



Franx % Ilhngworth 1996^ [Ziegler fc Bender 1997| ; [LonghetT^ 
et al. Iil98|; [Mehlert 199^ ; l0rgensen 1999[ ) . In contrast with 
high resolution index systems ( [Rose 1994| ; Jones fc Worthey 
199£j)7}hich promise a better separation of age and metalUi- 
tiy, the Lick/IDS system allows the investigation of dynam- 
ically hot galaxies that have intrinsically broad absorption 
lines. By plotting an age sensitive index and a metallicity 
sensitive index against each other, one can (partially) break 
the age /metallicity degeneracy and estimate, with the help 
of model predictions, the luminosity weighted age and metal- 
licity of an integrated stellar population (see Figure^). Most 
recently, j0rgensen (1999) used this methodology to inves- 
tigate the stellar populations of a large sample of early-type 
galaxies in the Coma cluster. She concluded that there are 
real variations in both the ages and the abundances while an 
anti-correlation between the mean ages and the mean abun- 
dances makes it possible to maintain a low scatter in scaling 
relations such as Mg-cro. CoUess et al. (1999) present similar 
conclusions from the analysis of a combination of the Mg- 
ctq relation and the Fundamental Plane in a large sample of 
cluster early- type galaxies. 

The spread in the ages for early-type galaxies and the 
anti-correlation of age and metallicity found by the previous 
authors supports the hierarchical picture for the construc- 
tion of galaxies in which galaxies form via several mergers in - 



volving star- formation (Baugh et al. 1996; Kauffmann 1996) 
However, the results are inconsistent with the conventional 
view that all luminous elliptical galaxies are old and coeval. 
In this picture the global spectrophotometric relations ob- 
served for ellipticals, for example the colour-magnitude re- 
lation (Visvanathan fc Sandagc 1977; Bower, Lucey fc Ellis 



1992; lerlevich 1998) are explained by the steady increase 
in the abundance of heavy elements with increasing galaxy 
mass. This increase arises naturally in galactic wind models 
such as that of Arimoto fc Yoshii (1987) and Kodama fc 
Arimoto (1997). 

Although with line-strength indices we can (partially) 
break the age/metallicity degeneracy this is by no means 
the last obstacle to overcome on our way to fully understand 
the stellar populations of early-type galaxies and the cause 
of scaling ralations. Since the late 70 's evidence has been 
accumulating that abundance ratios in galaxies are often 
non-solar. In particular the Magnesium to Iron ratio seems 
to be larger in luminous early -type galaxies when compared 
to solar neighbourhood stars ( O'Connell 197(:; Peletier 1989t 



Worthey, Faber fc Gonzalez 1992; Davics, Sadler fc Peletier 



1993 |; [Ifenry fc Worthey 1999^ [Jorgcnsen 1999| ). However; 
with only a very few exceptions (e.g., Weiss, Peletier fc Mat- 
teucci 1995), non-solar abundance ratios have not yet been 
incorporated in the model predictions. Among other issues 
this seems to be the most important single problem which 
prevents us from deriving accurate absolute age a nd metal- 
licit y est imates from integrated light spectroscopy ( Worthey 
199 ^). ^ Nevertheless with the current models and high S/N 
data we are able to study relative trends in ages and abun- 
dances as well as start to investigate the effe cts of non-solar 



In this paper, high S/N nuclear spectra of a complete 
sample of early-type galaxies in the Fornax cluster brighter 
than Mb ~ —17 are analyzed in the Lick/IDS system. The 
early results of this stu dy have already been presen ted in 
a letter to this journal (Kuntschner fc Davies 1998). This 
paper is organized as follows. Section ^ describes the sample 
selection and basic data reduction. The calibration of the 
line-strength indices to the Lick/IDS system is presented in 
Section ^ and the Appendix. Section ^ presents a consis- 
tency test for the model predictions and our measured line- 
strength indices. In Section ^ several index combinations are 
compared to model predictions. In particular the effects of 
non-solar abundance ratios, composite stellar populations 
and age/metallicity estimates of the integrated light are dis- 
cussed. In Section ^ observed index-uo relations are pre- 
sented. Relations between derived parameters such as age, 
metallicity and [Mg/Fe] with the central velocity dispersion 
are investigated in Section |^ We then discuss the implica- 
tions of our results in Section ^ and present the conclusions 
in Section ^. The fully corrected Lick/IDS indices of our 
sample are tabulated in the Appendix. 



2 THE OBSERVATIONS AND DATA 
REDUCTION 

2.1 The sample 

Our sample of 22 early-type galaxies has been selected 
from the catalogue of Forna:x galaxies (Ferguson 1989, here- 
after F89), in order to obtain a complete sample down to 
Bt = 14.2 or Mb = -17.[| We have adopted the morpho- 
logical classifications given by F89 and checked them with 
images we obtained on the Siding Spring 40" telescope. From 
these we noted a central dust lane in ESO359-G02 and a 
central disc in ES0358-G59 which led us to classify them 
as lenticular galaxies. NGC 1428 was not observed because 
of a bright star close to its centre. We also added the ellip- 
tical galaxy IC2006 to our sample, as it was not classified 
by F89. The bona-fide elliptical NGC3379 was observed as 
a calibration galaxy. The observations were carried out with 
the AAT (3.9m) on the nights of 1996 December 6-8 using 
the RGO spectrograph. The characteristics of the detector 
and the instrument set-up is given in Table |l]. 

Typically, exposure times were between 300 and 1800 
sec per galaxy (see Table ^ for a detailed listing). For most 
of the observations the slit was centred on the nucleus at 
PA = 90°. The seeing was generally better than one arc- 
sec. Additionally we observed 15 different standard stars 
(mainly K-giants) during twilight to act as templates for 
velocity dispersion measurements as well as t o calibrate our 



line-s trength indices to the Lick/IDS system (Worthey et al 



abundance ra tios for individual elements (Worthey 1998 
Peletier 2o6o|) . 



1994). The spectrophotometric standard stars GD 108 and 
L745-46A were observed to enable us to correct the contin- 
uum shape of our spectra. Table |^ lists all observed standard 
stars with their spectral types (obtained from SIMBAD, oper- 
ated by CDS, Strasbourg) and also comments on their use as 



t Adopting a distance modulus oim — M = 31.2: based on I-ban d 
surface brightness fluctuations ( Jensen, Tonry fc Luppino 199 j ). 
This corresponds to Ho ~ 80 km s~^ Mpc"^ for a flat Universe. 
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Table 1. The instrumental set-up 



Table 3. Log of observations; stars 



Telescope 

Dates 

Instrument 



AAT {3.9m) 

6-8 December 1996 

RGO spectrograph 



Spectral range 

Grating 

Dispersion 

Resolution (FWHM) 

Spatial Scale 

Slit Width 

Detector 

Gain 

Read-out- noise 
Seeing 



4243 - 5828 A 
600 V 

1.55 Apixel"^ 
~ 4.1 A 
0'.'77 pixel -1 
2'.' 3 

Teklk #2 (24 fim^ pixels) 
1.36 e-ADU"^ 
3.6 e~ rms 

~1" 



Table 2. Log of observations: galaxies 





Tvoe 


Bt 


Exp. time 


PA 






[ma,E?l 


[sec] 


L J 


NGC 1316 


SO pec 


9.4 


300 


90 


NGC 1336 


E4 


13.1 


1200 


90 


NGC 1339 


E4 


12.5 


600 


90 


MOO 10C1 

IN GO 1351 


E5 


12.5 


900 


90 


NGC 1373 


E3 


14.1 


1200 


90 


NGC 1374 


EO 


12.0 


900 


90 


NGC 1375 


SO 


13.2 


900 


90 


NGC 1379 


EO 


11.8 


600 


90 








600 


90 


NGC 1380 


SO 


10.9 


300 


90 


NGC 1380A 


SO 


13.3 


900 


90 


NGC 1381 


SO 


12.4 


480 


90 








900 


140 








900 


50 


NGC 1381" 


SO 




1800 


50 


NGC 1399 


EO, cD 


10.6 


300 


90 








600 


181 


NGC 1404 


E2 


11.0 


300 


90 








300 


90 


NGC 1419 


EO 


13.5 


1200 


90 


NGC 1427 


E4 


11.8 


600 


90 








1800 


79 








1800 


169 


IC 1963 


SO 


12.9 


600 


90 


IC 2006 


E 


12.2 


600 


90 


ESQ 359-G02 


SO 


14.2 


1200 


90 


ESQ 358-G06 


SO 


13.9 


1800 


90 


ESQ 358-G25 


SO pec 


13.8 


1200 


90 


ESQ 358-G50 


SO 


13.9 


1200 


90 


ESQ 358-G59 


SO 


14.0 


1200 


90 


NGC 3379^* 


El 


10.2 


300 


90 



Notes - "offset from nucleus; 'calibration galaxy, non Fornax 
member 

Lick/IDS standard, velocity standard or spectrophotometric 
standard. 



2.2 Basic data reduction 

Most of the basic data reduction steps liave been performed 
witii packages under IRAF. For eacli niglit individually tiie 
science frames were overscan corrected and a bias frame was 



Ncime 


lype 


cominent 


HD004656 


K4IIIb 


Lick/IDS std 


HD037160 


KOIIIb 


Lick/IDS std 


HD040657 


K1.5III 


velocity std 


HD047205 


Kim 


Lick/IDS std 


HD050778 


K4III 


Lick/IDS std 


HD054810 


KOIII 


Lick/IDS std 


HD058972 


K3III 


Lick/IDS std 


HD061935 


G9III 


Lick/IDS std 


HD066141 


K2III 


Lick/IDS std 


HD071597 


K2III 


velocity std 


HD083618 


K2.5III 


Lick/IDS std 


HD088284 


KOIII 


Lick/IDS std 


HD095272 


Kim 


Lick/IDS std 


HD219449 


KOIII 


Lick/IDS std 


HD221148 


K3III variable 


Lick/IDS std 


L745-46A 


DF 


SDOC. std flokc 1971) 


GD 108 


sd:B 


spec, std (|okc 199C ) 



subtracted. A few bad columns were removed by linear inter- 
polation. From several domeflats and skyflats a final flatfield 
accounting for tlie pixel-to-pixel variations and vignetting 
was constructed and applied to tiie frames. Cosmic rays were 
r emoved using the cleanest task in the REDUCEHE package 
( [Cardiel et al. 1998 ) . This task automatically detects and re- 



moves cosmic rays via a sophisticated deviation algorithm, 
while at the same time one can interactively inspect poten- 
tial cosmic rays in sensitive areas such as close to the galaxy 
centre. The wavelength solution was determined from Th- 
Ar-lamp spectra which were taken before and after most of 
the science observations. The rms residual in the wavelength 
fit was typically 0.1 - 0.2 A. Finally the sky was subtracted. 

The central spectrum for each galaxy was extracted by 
fitting a low order polynomial to the position of the centre 
along the wavelength direction, re-sampling the data in the 
spatial direction and finally co-adding the spectra within a 
5 pixels aperture yielding an effective aperture of 2'.'3x3'.'85. 
Multiple exposures of the same galaxy were combined. The 
resulting S/N in the spectra ranges from ~ 30 [A~^] for the 
faintest galaxies to more than 100 [A~^] for the brightest 
ones (measured in a ~100 A wide region just bluewards of 
the Mg h feature) . For stars we used the IRAF task apall to 
extract Id-spectra. All galaxy and stellar spectra were log- 
arithmically rebinned to a common wavelength range and 
increment. Finally the continuum shape of our spectra was 
corrected to a relative flux scale with the help of the spec- 
trophotometric standard stars. 

2.3 Kinematics 

In order to correct the line-strength indices for velocity dis- 
persion broadening and to construct index-cro relations we 
need to measure the central velocity dispersion for each 
galaxy. Estimates were derived with the Fourier correlation 
quotient (FCQ, version 8) method (Bender 1990; Bendi 
Saglia & Gerhard 1994). For the FCQ analysis the spectra 
were rebinned to twice the original spectral sampling and 
a wavelength range of 4876 to 5653 A was extracted. Note 
that the H/3 feature is excluded from the wavelength range 
as it proved to be a source of severe template mismatch for 
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galaxies with strong Balmer absorption. As we only consider 
central spectra in this paper we fit a pure Gaussian profile 
to the broadening function, neglecting higher order terms. 
To check the reliability of the FCQ analysis we used eight 
different G & K-giant template stars. For galaxies with a 
central velocity dispersion of ao > 70 km s~^, all eight tem- 
plate stars give very similar results and an average value was 
adopted. The rms scatter between different template stars 
is 0.007 in log units for galaxies with ctq > 100 km s~^. For 
galaxies with 70 < ao < 100 km s"'^ the rms scatter in- 
creases to 0.024 and for galaxies with ao < 70 km s~^ we 
find an rms scatter of 0.074. The uncertainty introduced by 
different template stars was comparable or larger than the 
internal error estimates of the FCQ program. Note, that for 
galaxies with ao < 70 km s^^ some template stars gave a 
poor fit to the broadening function and were excluded from 
the template sample. Only remaining measurements were 
averaged. 

Using this procedure, velocity dispersions as low as 
~ 50 km s~^ could be recovered, although systematic er- 
rors will start to dominate for ao < 90 km s~^. As our 
spectral resolution is rather low compared to velocity dis- 
persions of ~ 50 — 60 km s~^ we emphasize that for these 
faint galaxies our velocity dispersions are only rough esti- 
mates. The final velocity dispersion errors for galaxies with 
a > 70 km s"^ (Alogao = 0.022) were derived by a liter- 
ature comparison (see Appendix, Figure Al). For galaxies 
with ao < 70 km s~^ we adopt the mean rms scatter of the 
template stars (Alogao = 0.074). 



3 LICK/IDS CALIBRATION 

The wavelength range of our spectra covers 16 different line- 
strength indices, such as Mg2, H/3 and H7A, in the Lick/IDS 
system which is described in detail in Worthey (1994, here- 
after W94), Worthey & Ottaviani (1997, hereafter W097) 
and Trager et al. (1998). In the following analysis we use an 
updated version of the W94 models which is available from 
Dr. G. Worthey's home page. The updates affect only mod- 
els where [Fe/H] < —1.0. and are most noticable for the H/3 
index. For a recent study of the behaviour of the Balmer in- 
dices at low metallicities see Poggianti & Barbaro (1997) and 
Maraston, Greggio & Thomas (1999). Before one can com- 
pare the measured indices with model predictions by e.g., 
W94 and Vazdekis et al. (1996, hereafter V96), the measure- 
ments have to be carefully calibrated to the Lick/IDS sys- 
tem. Generally there are three effects to account for: (a) the 
difference in the spectral resolution between the Lick/IDS 
system and our set-up, (b) the internal velocity broaden- 
ing of the observed galaxies and (c) small systematic offsets 
caused by e.g., continuum shape differences. 

(a) In order to account for diffrences in spectral resolu- 
tion we broadened the spectra with a Gaussian of wave- 
length dependent width, such that the Lick/IDS resolu- 
tion was best matched at each wavelength (see Figure 7 
in W097). After this step our spectra should resemble very 
well the general properties of the original spectra obtained 
by the Lick group. 

(b) In a second step we need to correct the indices for 
velocity dispersion broadening. The observed spectrum of a 
galaxy is the convolution of the integrated spectrum of its 



stellar population(s) by the instrumental broadening and 
the distribution of line-of-sight velocities of the stars. These 
effects broaden the spectral features, in general reducing the 
observed line-strength compared to the intrinsic values. In 
order to compare the raw index measurements for galaxies 
with model predictions we calibrate the indices to zero ve- 
locity dispersion. Spectra of 15 different G9-K4 giant stars 
were first broadened to the Lick/IDS resolution and then 
further broadened using a Gaussian to a = 20 — 360 km s~^ 
in steps of 20 km . The indices are then measured for 
each star and a-bin and a correction factor, C(a), such that 
C(a) =lndex(0) /Index(a) is determined. 

Figure BS in the Appendix shows the dependence of 
the correction factor on a for all 16 indices. Note that for 
the molecular indices Mgi and Mg2 and the index H7f|T] the 
correction factor is defined as C(a) =lndex(0) - Index(a). 
The scatter in C(a) at 360 km s~^ was < 5% for all indices 
but H/3. 

It is worth looking in detail why the H/3 velocity dis- 
persion correction seems to be so insecure. The derived cor- 
rection factors are only useful if the stars used for the simu- 
lations resemble the galaxy spectra. In principle one might 
expect a dependence of the correction factor on line-strength 
- but most indices do not show such a behaviour. In fact H/3 
is the only index where we find a significant influence of line- 
strength on the correction factor at a given a. It turns out 
that stars which exhibit H/3-absorption < 1.1 A lead to cor- 
rection factors of C(a) < 1.0 and stars with H/3-absorption 
> 1.1 A imply positive corrections. In the Fornax sample 
there are no galaxies with H/3 absorption line-strength of 
less than 1.4 A, hence only stars with a H/3 index greater 
1.1 A have been used to evaluate the correction factor (scat- 
ter < 5%). 

Another way to check the accuracy of the velocity dis- 
persion corrections is to use galaxy spectra with small in- 
ternal velocity dispersions as templates and treat them in 
the same way as stars. The galaxies NGC1373, NGC1380A, 
NGC1336, IC1963 and ES0358-G59 were used for this pur- 
pose. They span a range in H/3 absorption of ~ 1.7 — 3 A and 
a range in central velocity dispersion of ao — 54 — 96 km s~^. 
In Figure B2 the galaxies are represented by open circles and 
they agree very well with the stellar correction for most of 
the indices. As expected for H/3, the galaxies match the re- 
sults from stars with a H/3 absorption > 1.1 A. 

The final correction factors are derived by taking the 
mean of 15 stars and the five galaxies in each a-bin (solid 
line in Figure B2). The velocity dispersion corrections are 
applied by a FORTRAN program which reads in the raw index- 
measurements from continuum corrected and resolution cor- 
rected galaxy spectra. For each galaxy and index it applies a 
correction for velocity dispersion. The program linearly in- 
terpolates between a-bins and also adds the error from the 
velocity-dispersion correction factor to the raw Poisson er- 
ror of the spectra. As the error in the correction factor is 
much bigger than any error caused by uncertainties in a, we 
assumed the velocity dispersion of the galaxies to be error 
free. 

(c) Although we have matched very well the spectral 



t This index is actually not a molecular index but typical index 
values are close to zero hence a correction factor can degenerate. 
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Table 4. Lick/IDS Offsets 



Index 


offset (Lick/IDS - AAT) 


G4300 


+0.21 ±0.09 A 


Fe4383 


+0.60 ±0.13 A 


Ca4455 


+0.37 + 0.06 A 


Fe4531 


+0.00 + 0.10 A 


C24668 


-0.19 + 0.17 A 


H/3 


-0.05 + 0.04 A 


Fe5015 


+0.00 + 0.08 A 


A/To--, 




Mg2 


+0.023 + 0.003 mag 


Mgfe 


+0.15 + 0.09 A 


Fe5270 


+0.07 + 0.05 A 


Fe5335 


+0.00 + 0.08 A 


Fe5406 


+0.00 + 0.04 A 


Fe5709 


+0.00 + 0.06 A 


H7A 


+0.45 + 0.28 A 


H7P 


+0.00 + 0.14 A 



after the correction for velocity dispersion. Note that the 
Lick/IDS-system offset-error is a constant value and does 
not depend on the velocity dispersion of the galaxy itself. 
Therefore we did not include this error in the individual 
index errors but rather quote for each index a common 
offset error (see also Table The final corrected central 
(2'.'x3'.'85) index measurements and associated errors for the 
Fornax galaxies and NGC3379 are presented in Table D2 in 
the Appendix. For each galaxy we give the index measure- 
ment in the first row and the lo error in the second row. 

Note that for the galaxies NGC1381 and NGC1427 we 
combined three exposures yielding a very high S/N spec- 
trum. Here our index-error estimation taking into account 
only the Poisson error becomes invalid because of other er- 
ror sources such as the wavelength calibration, contirmum 
correction and aperture effects. By comparing individual ex- 
posures we established that 1.5 x the original Poisson error 
estimate is a good indicator of the true error. This adjusted 
error was adopted in Table D2 and for any further analysis. 



resolution of the Lick system, small systematic offsets of the 
indices introduced by continuum shape differences are gen- 
erally present (note that the original Lick/IDS spectra are 
not flux calibrated). To establish these offsets we compared 
our measurements for stars in common with the Lick/IDS 
stellar library. In total we observed 13 different Lick/IDS 
stars. Figure Bl in the Appendix shows the difference be- 
tween Lick/IDS measurements and ours after the mean off- 
set has been removed. The mean offsets and associated er- 
rors for each index are summarized in Table ^. The star 
HD221148 was excluded from the offset analysis because 
our index measurements proved to be very different from 
the original Lick/IDS measurements - possibly due to its 
variable nature (see Table |^ . The formal error in the offset 
is evaluated by the mean standard deviation of stars with 
respect to the mean offset divided by y^Ustars — 1. 

Most off the indices show small offsets to the Lick/IDS 
system, similar to the ones quoted in ( W097, Table 9) . The 
rather large offset in Mg2 is due to a well known difference 
in continuum shape. 

Recently Trager et al. (1998) published the Lick/IDS 
library of extragalactic objects including 7 galaxies in the 
Forna:x cluster and NGC3379. We can check our previous 
offset evaluation by comparing our galaxy measurements 
with Trager et al. For this purpose we extracted a 3 pixels 
central aperture (2'.'3x2'.'44) for our galaxies matching the 
Lick/IDS standard aperture of l'.'4x4". Our indices are then 
corrected for velocity dispersion as described in paragraph 
(b) and the offsets from Table |^ are applied. The results are 
overplotted in Figure 



Bl 



in the Appendix (filled symbols) 
The galaxies show for all indices more scatter around the 
mean offset than the stars which is somewhat reffected in 
the bigger error bars, but there are also some outliers. This 
is not surprising as seeing effects and aperture differences 
will introduce some non-reproduceable offsets for individual 
galaxies. Furthermore we note that the Lick group had to 
observe the Fornax galaxies at a very high airmass. With the 
possible exception of the indices G4300 and Fe4383 the off- 
set inferred from the galaxy comparison are consistent with 
the stellar comparison. 



4 A CONSISTENCY TEST OF THE MODEL 
PREDICTIONS 

For the following analysis of the nuclear stellar populations 
(Section ^) it is extremely important that our index mea- 
surements are accurately calibrated onto the Lick /IDS sys- 
tem which i s based on the Lick/IDS stellar library (Worthey 
et al. 1994). Here we investigate the accuracy and consis- 



tency of our calibration and the model predictions by pre- 
senting index-index plots which are almost degenerate in 
age and metallicity. In that way the model predictions cover 
only a small "band" of the parameter space and they should 
trace the relation of the galaxies if the models describe accu- 
rately the galaxy properties and our calibration is accurate. 
Figure ^ shows the relation between <Fe>^, Fe5015 and 
Fe5406 for our sample of galaxies (filled circles) and the orig- 
inal Lick/IDS sample of galaxies (Trager et al. 1998, small 
dots). Overplotted are models by W94 (black lines) and V96 
(grey lines). The plots show a good agreement between in- 
dex measurements and the model predictions. The reduced 
Poisson noise of our data set compared to the Lick/IDS mea- 
surements can be clearly seen (see also Figure caption). We 
note that the model predictions of Vazdekis and Worthey 
are in good agreement. 

A similar analysis of the three Mg indices is shown in 
Figure ^. Here we find a significant deviation of the measured 
index values compared to the model predictions for metal 
rich and/or old stellar populations. The deviations are seen 
in the Fornax sample as well as in the original Lick/IDS 
galaxy spectra (see also Worthey 1992, Figure 5.12 & 5.13). 
We therefore note that this discrepancy is inherent to the 
Lick/IDS system & models and any models which use the 
Lick/IDS fitting functions are likely to show the same offset. 
In Figure^ we present the Balmer indices H/3, H7A and H7F. 
Here we find generally good agreement with small deviations 
between model predictions and data at low values of H7A vs 
H7F which is present in our data and the original Lick/IDS 
measurements. 



The offsets listed in Table ^ were applied to all indices § <Fe>=(Fe5270+Fe5335)/2 
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Figure 1. Index vs index plots for three well established Fc- 
indicators in the Lick/IDS system. The filled circles and small 
dots represent AAT and Lick/IDS galaxy measurements respec- 
tively. The error bar in the upper left corner represents the av- 
erage observational error for Lick/IDS galaxies whereas for the 
AAT data the observational errors are shown for each individ- 
ual galaxy. The error bar in the lower right corner shows the 
rms error in the offset to the Lick/IDS system for the AAT data. 
Overplotted are model predictions by Worthey (1994, black lines) 
and Vazdekis et al. (1996, grey lines). Note that Worthey models 
use a Salpeter IMF whereas Vazdekis models use a bimodal IMF 
which is very similar to Salpeter for M > 0.6 Mq. 
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Figure 2. Index vs index plots for the three Mg indices in the 
Lick/IDS system. Overplotted are model predictions by Worthey 
(1994, black lines) and Vazdekis et al. (1996, grey lines). 
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Figure 3. Index vs index plots for three Balmer line indices in the 
Lick/IDS system. Symbol definitions are the same as in Figure ^ 
except for the galaxy ES0358-G25 which is represented by an 
open triangle. This galaxy is affected by emission in the Balmer 
lines. Overplotted are model predictions by Worthey (1994, black 
lines), Worthey & Ottaviani (1997, black lines) and Vazdekis et al. 
(1996, grey lines). 

Figures § to ^ suggest that our Lick/IDS calibration is 
very consistent with the original galaxy measurements of the 
Lick group. However, we note that small, systematic offsets 
exist between the parameter space covered by galaxies and 
the model predictions for Magnesium at high index values 
and for Balmer lines at low index values. 



5 THE NUCLEAR STELLAR POPULATIONS 

The aim of this section is to derive estimates of the mean 
(luminosity weighted) ages and metal abundances of early- 
type galaxies in the Fornax cluster. As pointed out by W94, 
the determination of the ages and metallicities of old stellar 
populations is complicated by the similar effects that age 
and metallicity have on the integrated spectral energy dis- 
tributions. However, this degeneracy can be partially broken 
by plotting a particular age sensitive index, such as one of 
the Balmer line indices, against a more metallicity sensi- 
tive index. The usefulness of this approach has been demon- 
strated by many authors ( Glonzalez 1993 ; [Fisher, Franx & 
[Uingworth 1995| ; [Mehlert 1998t [Kuntschner fc Davies 199^ ; 
j0rgensen 199S;). However, as we will see in this Section, 
among other issues the treatment of non-solar abundance 
ratios is a crucial parameter in the determination of abso- 
lute age and metallicity estimates. We will also investigate 
the effects of nebular emission and composite st ella r p op- 



5.3 



ulations on age/metallicity estimates in Sections 5.2 & 
respectively, before we present our best age / metallicity es- 
timates of the Fornax early-type galaxies in Section 5.4 
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Figure 4. The two age sensitive indices H/3 and H7A plotted against six metallicity indicators: Mg2, C24668, Ca4455, Fe3, <Fe> and 
Fe5406. Models by Worthey (1994, black lines), Worthey & Ottaviani (1997, black lines) and Vazdekis et al. (1996, grey lines) are 
overplotted. The solid lines represent isoage lines whereas the dashed lines are lines of constant metallicity. Filled circles and open circles 
represent ellipticals and SOs respectively. The star and open triangle represent possible post-starburst and starburst galaxies respectively. 
The arrow attached to ES0358-G25 (open triangle) indicates an emission correction. The symbol size is scaled with the central velocity 
dispersion of the galaxies. The cross in the upper right corner of each panel indicates the rms uncertainty in the transformation to the 
Lick/IDS system. Observational errors are plotted as individual error bars on the data points. 
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5.1 Non-solar abundance ratios 

In Figure ^ we present age/metallicity diagnostic diagrams 
of six metallicity sensitive indices (Mg2, C24668, Ca4455, 
Fe3, <Fe> & Fe5406) plotted against the age sensitive 
Balmer line indices H/3 and Hta (the new index Fe3 is 
defined in Equation |l]). Figure ^ is a reproduction from 
Kuntschner & Davies (1998) with minor data up-dates. 
Overplotted are model predictions from W94, W097 (black 
lines) and V96 (grey lines) . Solid lines represent lines of con- 
stant age and the dashed lines are lines of constant metal- 
licity. The Worthey models span a range in age of 1.5-5 
Gyr with [Fe/H]=-0.225 to 0.5 and 8-17 Gyr with [Fe/H]=- 
2 to 0.5. The V96 models span a range in age of 1-17.4 Gyr 
with [Fe/H]=-0.7 to 0.4. The direction of increasing age and 
metallicity is indicated in Figure ^a by arrows. 

Our previous result from Kuntschner & Davies (1998), 
that Fornax ellipticals form a sequence in metallicity at old 
ages and that the SOs spread to younger ages is confirmed 
in all diagrams. Examining the diagrams in detail one can 
see that the mean age and metallicity of the sample changes 
from diagram to diagram; e.g., the ellipticals appear older 
and more metal poor in the <Fe> vs H7A diagrams com- 
pared to the Mg2 vs H7A diagram. This effect was previ- 
ously reported and recently reviewed by Worthey (1998). It 
is now widely accepted that this discrepancy in the model 
predictions is caused by non-solar abundance ratio effects. 
For example Mg as measured by the Mg2 index is over- 
abundant co mpared to Fe in luminous ellipt i cal galaxies, i.e., 
[Mg/Fe ] > (lO'Connell 1976|; [Peletier 1989|; [Worthcy, Fab~ 
& Gon 
Peletiei 



alez 1992|; [Davies, Sadler fc Peletier 1993t [Weiss 



fc Matteucci 1995| ; [l0rgensen 1997| ; p0rgcnsen 1999D 



The Mg over abundance can be examined in a M g-index 
vs Fe- index plot ( Worthey, Faber fc Gonzalez 1992 ). In such 
a diagram the model predictions cover only a narrow band 
in the parameter space as effects of age and metallicity are 
degenerate. Figure M shows plots of <Fe> & Fe5270 vs Mg2 
for the Fornax sample. Overplotted are model predictions 
from W94, V96 and Weiss et al. (1995). We assume that the 
models refiect solar abundance ratios if not stated other- 
wise, i.e., [Mg/Fe] — 0. If the model predictions accurately 
resemble the galaxy properties they should trace the ob- 
served relation. The measured line-strength of most of the 



SOs agrses with the model predictions, perhaps 3-4 galax- 
ies having slightly low Mg2 absorption compared to Fe5270 
& <Fe>. However, for most of the ellipticals and the SO 
NGC138Q, the models predict too little Mg-absorption at a 
given Fe-absorption strength. Additionally the most metal 
rich galaxies are the furthest away from the model grids. Us- 
ing the Mg-overabundance correction by Greggio (1997, see 
Figure ^a) and the models for [Mg/Fe]=0.45 by Weiss et al. 
(1995), we conclude that the stellar populations of Fornax 
ellipticals and the bulge of NGC1380 are Mg-overabundant 
relative compared to Fe. The overabundance ranges between 
[Mg/Fe]=0.0 to ~0.4. We note that there exists considerable 
spread in overabundance at a given Fe-line strength in our 
sample. 

Non-solar abundance ratios do exist not only in ellipti- 
cal galaxies but also in our own galaxy where stars show an 




Figure 5. (a) Mg2 vs Fe5270 equivalent width diagram for the 
complete sample of Fornax early-type galaxies. Overplotted are 
models by Worthey (1994) and a correction for [Mg/Fe]=0.4 for 
the 17 Gyrs isoage line (taken from Greggio 1997). (b) Mg2 vs 
<Fe> diagram. Overplotted are models by Vazdekis (1996, grey 
lines) and two models by Weiss, Peletier & Matteucci (1995, dot- 
dashed lines). The Weiss et al. models are calculated for three 
ages, 12, 15 & 18 Gyrs (dot-dashed lines represent lines of con- 
stant age) at Z=0.02, 0.04 & 0.07, a mixing length parameter 
Q^MLT = 1-5 and somewhat different mixes of heavy elements. 
Steps in metallicity are shown as diamonds for [Mg/Fe] = 0.0 
and as triangles for [Mg/Fe] = 0.45. Symbol definitions as in Fig- 



overabundance for a-elements|^ at [Fe/H]^0.0 (Edvardsson 



et al. 1993 



; |McWilliam 199^ ). Of course if those stars are in- 
corporated in a stellar library which in turn is used for model 
predictions of integrated stellar populations, the predictions 
will be somewhat a-element overabundant at low metallici- 
ties. We therefore note that models which use the Lick/IDS 
fitting functions are probably a-element overabundant at 
low metallicities which makes it more difficult to interpret 
trends in diagrams such as Figure 

Several indices covered by our wavelength range show 
deviations from the model predictions when compared to 
the average Fe index: Mgi, Mg2, Mg &, Fe5709 & C24668. 
Fe5709 is a very weak index and its correction for velocity 
dispersion broadening may well be insecure, so we cannot 
draw any firm conclusions. C24668 is an important index 



a includes the elements O, Mg, Si, S, Ca, and Ti 
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because it shows the stro ngest total met aUicity sensitivity 
in the Lick/IDS system (Worthey 1998) and is therefore 



preferentially used in age/metallicity diagnostic diagrams. 
In Figure ^ we present a plot of C24668 vs Fe3. Fe3 is a 
combination of three prominent Fe lines, thus maximizing 
its sensitivity to Fe while minimizing the Poisson errors: 



Fe3 = 



Fe4383 + Fe5270 + Fe5335 



(1) 



As a consequence of the extreme metaUicity sensitiv- 
ity of C24668, the models are not as degenerate as in the 
previous plots. Nevertheless it is clear that for a C24668 ab- 
sorption strength in excess of ~6 A the model predictions 
do not follow the observed trend (see also Kuntschner 1998) . 
Hence, we conclude that C24668, or better at least one of the 
species contributing to the index, is overabundant compared 
to Fe in metal rich Fornax galaxies. Can this overabundance 
be caused by Mg as seen in the Fe vs Mg plot (Figure ^)? 
Due to the proximity of metal absorption lines in the opti- 
cal wavelength region non of the Lick/IDS indices measures 
the abundance of only a particular element such as Fe or 
Mg. There are always c ontributions from othe r elements or 
molecules to an index (Tripicco & Bell 1995). In particu- 



lar, the C24668 index has a relatively wide central bandpass 
(86.25 A) including a wide range of metal lines. The most 
dominant species here is carbor i in form of C2-bands which 
blanket the central bandpass ( Tripicco fc Bell 1995 ). Yet, 
more important here is that according to Tripicco & Bell 
(1995; Table 6, cool giants) the C24668 index decreases when 
the Mg-abundance (or Oxygen-abundance) is increased (at 
fixed abundances of all other elements) . Therefore the 'over- 
abundance' of C24668 cannot be caused by Mg. What ex- 
actly drives the overabundance of C24668 compared to Fe 
remains to be seen. 

The overabundance of certain elements compared to Fe 
has a profound effect on the use of age/metallicity diag- 
nostic diagrams if the model predictions reflect only solar 
abundance ratios. Recalling Figure ^ (see also Kuntschner 
& Davies 1998) we find that not only are the metallici- 
ties (measured as [Fe/H]) overestimated due to the C24668 
overabundance but furthermore much of the age upturn at 
high metallicities is caused by the overabundance. The ef- 
fect of changing the age estimates is caused by the residual 
age/metallicity degeneracy which is still present in all dia- 
grams in Figure ^. Only if the index-combination breaks the 
degeneracy completely, i.e., if lines of constant age and con- 
stant metaUicity are perpendicular, would the age estimates 
not be affected by non-solar abundance ratios. We further 
note that trends in abundance ratios within a dataset such 
as our Fornax sample (increasing Mg/Fe with galaxy mass) 
can lead to artificial relative age trends in diagrams such 
as in Figure Taking into account not only the uncertain- 
ties introduced by non-solar abundance ratios but also other 
model parameters such as which isochrone library to use, it 
seems very insecure to derive absolute age estimates from 
the currently available stellar population models. 

Introducing non-solar abundance ratios in model pre- 
dictions is rather complicated, as accurate model predic- 
tions do need a stellar library covering the whole parameter 
space of Te, log g, [Fe/H] and [Mg/Fe]. Furthermore, new 
isochrone calculations may well be needed for each [Mg/Fe] 
bin: recently Salaris & Weiss (1998) suggested that scaled- 




4 6 
C,4668 |A| 



Figure 6. C24668 equivalent width vs Fe3 equivalent width. 
Overplotted are models by Worthey (1994). The C24668 index 
shows evidence of overabundance compared to Fe3 at strong ab- 
sorption strength (> 6.5 A). The symbol size is scaled with the 
central velocity dispersion of the galaxies. 



solar isochrones cannot be used to replace Mg-enhanced ones 
at the same total metaUicity. The latter will not only change 
the model predictions for indices such as Mg2 but may affect 
all indices and in particular the age sensitive ones such as 
H/3 and H7A (see also Worthey 1998). However, note that 
Weiss, Peletier & Matteucci (1995) concluded in their study 
that scaled solar isochrones are sufficient to calculate model 
predictions for non-solar abundance ratios. 

Another way to examine non-solar abundance ratios is 
to compare the metaUicity estimates derived from different 
metal lines using the same age indicator. Figure |^ com- 
pares the metaUicity estimates taken from Mg2, C24668, 
Fe5406, Ca4455 vs H/3 diagrams with the estimates taken 
from a Fe3 vs H/9 diagram. Here Fe3 serves as our mean Fe- 
abundance indicator. The metaUicity estimates are derived 

from the V96 model^ To get more accurate estimates, the 
age/metallicity-grid was expanded to a step size of 0.025 
in [Fe/H] by linear interpolation. Furthermore the diagram 
was extrapolated to [Fe/H]=-|-0.7 by linear extrapolation. 
The age range of 1 to 17.4 Gyrs is covered by 18 grid points. 
Errors on the metaUicity estimates were derived by adding 
and subtracting the index error for each galaxy individually 
(Poisson error and Lick/IDS offset error added in quadra- 
ture) and re-deriving the metaUicity estimates. The final un- 
certainty displayed in Figure ^ was taken to be 0.7 times the 
maximum change in [Fe/H]. 



II These models have a bimodal IMF which is very similar to 
Salpeter for M > 0.6 Mq. Note that for an age of ~17 Gyrs V96 
models predict 0.1 — 0.2 A less H/3-absorption compared to W94 
models. This of course will affect the absolute age estimates but 
has little affect on the metaUicity estimates. 
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Figure 7. Metallicity estimates derived from four age-metallicity diagnostic diagrams all using H/3 as age indicator but different metal 
lines (Mg2, Fe5406, C24668 & Ca4455) are compared with tfie metallicity estimates from the Fe3 vs H/3 diagram. The filled circles 
represent elliptical galaxies and the open circles stand for the SOs. The symbol size is scaled with the central velocity dispersion of the 
galaxies. 



In panel (a) of Figure |^ we can clearly see that for el- 
liptical galaxies Mg2 gives metallicity estimates which are 
larger than those derived from Fe3 and there is a trend that 
the Mg-overabundance increases with increasing metallicity. 
Most of the SOs are consistent with solar or slightly less than 
solar abundance ratios of Mg/Fe. However, the (more lumi- 
nous) SOs NGC1380 fc NGC1381 show a weak overabun- 



dance c f Mg. The index C24668 (panel b) gives on average 
high metallicity estimates compared to Fe3. Although three 



galaxies with just above solar metallicity show solar abun- 
dance ratios. As expected the Fe-index Fe5406 (panel c) is 
in good agreement with the estimates derived from FeS. The 
Ca4455 index (panel d) gives marginally higher metallicity 
estimates compared to the FeS indicator. We note that the 
Ca4455 index is more sensitive to a mix of hea vy elements 
than to C alcium on its own despite its name (Iripicco & 



BeU 1995) 



In conclusion we can confirm our previous results that 
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Mg and C24668 are overabundant compared to Fe. The Mg- 
overabundance follows a trend where metal rich (and lumi- 
nous) Es show a stronger overabundance than less luminous 
and metal poor galaxies. 

5.2 Nebular emission in early- type galaxies 

So far we have concentrated on breaking the age/metallicity 
degeneracy and the treatment of non-solar abundance ratios. 
A further important issue when estimating ages and metal- 
licities from line strength indices is nebular emission. Ellip- 
tical galaxies normally contain much less dust and ionized 
gas than spirals, in fact, for a long time they were regarded 
as dust and gas free. However, spectroscopic surveys of large 
samples of early-type galaxies revealed that about 50-60% 
of the galaxies show wea k optical emission lines ( Phillips 



et al. 1986; Caldwell 1984). Typically the reported strength 



of emission lines such as [Oil] , [Ha] and [Nil] A6584 indicates 
the presence of only 10^ — 10^ of warm ionized gas in the 
centre. A more recent study of 56 bright elliptical galaxies 
by Goudfrooij et al. (1994) detected ionized gas in 57% of 
their sample and confirms the amount of ionized gas present. 
Additionally, HST images of nearby bright early-type galax- 
ies revealed th at approximately 70-80% sh ow dust features 
in the nucleus ( van Dokkum fc Franx 1995| ). Stellar absorp- 



tion line-strength measurements can be severely affected if 
there is emission p resent in the galaxy which we akens the 
stellar absorption (Goudfrooij & Emsellem 1996). For ex- 
ample, nebular H/3-emission on top of the integrated stellar 
H/3-absorption weakens the H/3-index and leads therefore to 
wrong, i.e., too old age estimates. 

The spectrum of ES0358-G25 shows clear emission in 
H/3 and H7 along with weak [O III] emission (see Kuntschner 
& Davies 1998, Figure 3). As a consequence, the age is 
overestimated in Figures ^ (a)-(f). The arrow attached to 
ES0358-G25 indicates a rough emission correction. How- 
ever, it is extremely difficult to accurately correct the H/3- 
index in individual galaxies for emission contamination. A 
much better method to reduce emission contamination is to 
use higher order Balmer lines such as H7 as t hey are less 
affected by nebular emission ( lOstcrbrock 1989| ). Indeed, in 
Figures ^ (g)-(l) the galaxy ES0358-G25 moves to much 
younger ages. As none of the other galaxies move signifi- 
cantly to younger ages we conclude that nebular emission is 
not very prominent in our Fornax sample. This is supported 
by the absence of strong [OIII]A5007 emission. Only 5 galax- 
ies show emission above our detection limit of ~0.2 A. The 
strongest emission is detected in ES0358-G25 with 0.7 A 
equivalent width (for details see Kuntschner 1998). 

5.3 Effects of composite stellar populations 

Most of the SOs in our sample have luminosity weighted 
young stellar populations with some of them having also 
high metallicities when compared to single-burst stellar pop- 
ulation (SSP) models. However, these galaxies show only a 
central young stellar population on top of an underlying 
older one as opposed to be entirely young. It is not straight 
forward to compare composite ste llar populations with SSP 
models (ie Jong and Davies 1997). So, how reliable are the 
age, metallicity and abundance ratio estimates taken di- 
rectly from SSP models for these young galaxies? In order to 



explore this issue we calculated model predictions based on 
V96 for simple composite stellar populations. Two represen- 
tative tracks are shown in Figure H for two age/metallicity 
diagnostic diagrams and a plot of Fe3 vs Mg2 in order to 
explore the behaviour of abundance ratios: Model A is a 
15 Gyr old (90% in mass) stellar population plus a burst 
(10% in mass) of varying ages from 0.1 Gyr to 3 Gyrs. 
Both populations have solar metallicity. For Model B we re- 
duced the burst fraction to 1% (in mass) while the other 
parameters are the same as in Model A. 

Overall one can see in the age/metallicity diagnostic 
diagrams that for a short time the burst population will 
dominate the integrated light leading to strong H/3 absorp- 
tion and weak metal-line absorption. Then the underlying 
old population becomes more and more important and after 
~3 Gyrs the galaxy is almost back to its original place in the 
diagram. However, the burst strength influences the exact 
track which the galaxy takes in the diagram. For a burst of 
10% or 20% (not shown) in mass the tracks follow roughly 
the solar metallicity line in the normal SSP models. Yet, for 
a small burst (1% in mass) the integrated light looks for a 
short while as having metallicities well above solar. This ef- 
fect is more pronounced for Mg2 than for Fe3. Of course, this 
in turn leeds to an artificially created overabundance when 
these galaxies are compared to SSP models (see Figure |^). 
For bursts stronger than a few percent the abundance ratios 
are not significantly affected. 

In summary we find that composite stellar populations 
and in particular small (in mass) bursts, such as used in our 
simple models, can lead to an overestimation of the metallic- 
ity in the context of SSP models. Abundance ratios can be 
affected in the sense that the Mg/Fe ratio is too strong. Our 
model calculations show that these conclusions qualitatively 
hold if the metallicity is changed or different metallicities 
are combined. A more thorough investigation of these issues 
would be very valuable but is beyond the scope of this pa- 
per (see Hau, Carter & Balcells (1999) for a more detailed 
analysis) . 



5.4 Best age and metallicity estimates 

Having examined some of the fundamental problems with 
applying stellar population model predictions to observed 
line-strength indices we present in Figure ^ what we con- 
sider our best age/metallicity diagnostic diagram. A mean 
Fe-index (Fe3) is plotted against an emission robust higher 
order Balmer line (H7A). Due to the lack of model predic- 
tions with non-solar abundance ratios we decided to avoid 
indices which are affected by overabundance problems (e.g., 
Mg & C24668). Instead we use here a combination of Fe- 
indices (Fe3) as metal indicator which will bias our results 
towards the Fe abundance. We note however, that our metal- 
licities are not to be understood as total metallicity but 
rather as a good estimate of the Fe-abundance. Any non- 
solar abundance ratios which affect H7A are ignored. Model 
predictions by W94 and V96 are overplotted in Figure ^. 

The ellipticals form a sequence of metallicity at roughly 
constant age. The centres of the bright SOs NGC1380 & 
NGC1381 follow the sequence of Es. The remaining SOs cover 
a large range in metallicity and spread to much younger 
luminosity weighted ages than the Es. We emphasize that 
these age and metallicity estimates are central luminosity 
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Figure 8. Evolutionary tracks based on Vazdekis et al. (1996) are sliown in tliree index-index diagrams for a composite (Model A) of 
a 15 Gyrs old (solar metallicity, 90% mass) and a young stellar population (solar metallicity, 10% mass) at burst ages 0.1, 0.13, 0.16, 
0.2, 0.25, 0.32, 0.4, 0.5, 0.63, 0.79, 1.0, 1.26, 1.58, 2.0, 2.51 & 3.16 Gyrs. The plus symbols along the tracks indicate the time steps. Age 
steps of 0.32 and 1 Gyrs are also indicated by numbers in panel (a). Model B represents a burst of 1% (in mass) strength with the same 
metallicities and age steps as in Model A. Age steps of 0.32, 0.63 and 1.26 Gyrs are indicated in panel (a) & (c). The region of normal 
SSP models (Vazdekis et al. 1996) is shown as thin lines. 
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Figure 9. Fe3 equivalent width vs H-ja equivalent width. Filled 
circles and open circles represent ellipticals and SOs respectively. 
The star and open triangle represent possible post-starburst and 
starburst galaxies respectively. The cross in the upper right corner 
of each panel indicates the rms uncertainty in the transformation 
to the Lick/IDS system. The symbol size is scaled with the central 
velocity dispersion of the galajdes. Note that the two bright SOs 
are somewhat hidden in the sequence of Es. 



weighted estimates and for apparently young galaxies the 
derived parameters are somewhat more insecure (see previ- 
ous discussion about the effects of composite stellar popula- 
tions) . The age and metallicity gradients within the galaxies 



will be discussed in a future paper. Fornax A, a bright pecu- 
liar SO, shows strong Balmer lines and strong Fe absorption 
which translates into a luminosity weighted young and metal 
rich stellar population. As we will see in the next section all 
other young or metal poor SOs have velocity dispersions of 
o"o ^ 70 km s~^. The two galaxies with the weakest metal 
lines and strong H/3 & H7A absorption (ESO359-G02, cross 
and ES0358-G25, open triangle) appear to be different from 
the rest of the sample. These galaxies are likely to be post- 
starburst or starburst galaxies respectively. They have re- 
markable spectra for early-type galaxies, showing blue con- 
tinua, strong Balmer lines, and weak metal lines. These 
galaxies are amongst the faintest in our sample and are ~3° 
away from the centre of the cluster (see also Kuntschner & 
Davies 1998). 



LINE-STRENGTH INDICES AND THE 
CENTRAL VELOCITY DISPERSION 

The central velocity dispersion ao of earl y- type galaxies is 



known to correlate strongly with colours (Bower, Lucey & 



Ellis 1992) and the absorption strength of the Mg-absorption 
feature at 5174 A ([Terlevich et al. 198l|; |Burstein et al . 198^ 



CoUess 



Bende r, Burstein fc Faber 199^ ; j0rgensen 1997; 

1999| ). The relatively small scatter about these rela- 



et al. 



tions imply that the dynamical properties of galaxy cores are 
closely connected with their stellar populations. However, 
analysing the Mg-ao relation for a sample of 736 mostly 
early-type g alaxies in 84 clusters, the EFAR group (CoUess 



et al. 1999) finds a rather large dispersions in age (40%) 
and in metallicity (50%) at fixed velocity dispersion using 
the constraints from the Mg-cro relation and the Fundamen- 
tal Plane. Correlations of other metal indices, such as <Fe>, 
with the central velocity dispersion have long been expected 
but so far relations have shown a large scatter and only weak 
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Table 5. Scaling relations 



Table 6. Scaling relations - comparsion of slopes with j0rgensen 
(1997,1999) 



index 










scatter 
[mag] 


Mg2 


= 


(0.191 ±0.023) log 


(TO - 


(0.127 ±0.054) 


0.017 


Mgi 




(0.136 ±0.015) log 


(TO - 


(0.158 ±0.035) 


0.014 


Mgfe' 




(0.102 ±0.020) log 


(TO - 


(0.056 ± 0.044) 


0.011 


C24668 


= 


(0.090 ± 0.018) log 


(TO - 


(0.110 ±0.042) 


0.012 


Fe3' 


= 


(0.038 ± 0.011) log 


O-Q + 


(0.014 ±0.025) 


0.005 


Fe4383' 


= 


(0.043 ± 0.019) log 


O-Q + 


(0.037 ± 0.045) 


0.007 


Fe4531' 




(0.036 ± 0.010) log 


fTO + 


(0.009 ±0.023) 


0.007 


Fe5015' 




(0.036 ± 0.008) log (TO + 


(0.002 ±0.019) 


0.005 


Fe5270' 




(0.029 ± 0.009) log 


fTO + 


(0.024 ± 0.020) 


0.004 


Fe5335' 




(0.043 ± 0.009) log 


(TO - 


(0.017 ±0.020) 


0.005 


Fe5406' 




(0.023 ± 0.012) log 


(TO + 


(0.023 ±0.026) 


0.005 


Ca4455' 




(0.035 ± 0.017) log 


(TO + 


(0.014 ±0.038) 


0.009 


H/3' 




-(0.020 ±0.007) log 


CO + 


(0.106 ±0.015) 


0.004 


H7a' 




-(0.045 ±0.019) log 


(TO - 


(0.038 ± 0.044) 


0.010 


H7f' 




-(0.049 ±0.017) log 


O"0 + 


(0.018 ±0.037) 


0.009 



index 



this data 



literature 



reference 



corr elations ( Fisher, Franx fc Illingworth 199( ; j0rgensen 



1997; Jargensen 199!:). However, we will demonstrate that 
galaxies in the Fornax cluster do show a clear correlation 
between Fe-indices and central velocity dispersion. 

Following CoUess et al. (1999) we find it more conve- 
nient to express the "atomic" indices in magnitudes like the 
"molecular" index Mg2. The new index is denoted by the 
index name followed by a prime sign ['], e.g., Mg6'. Note 
that by using only the logarithm of the atomic index, one 
introduces a non linear term in comparison to the magni- 
tude definition. Furthermore negative index values such as 
for the H7A index cannot be put on a simple logarithmic 
scale. A priori it is not clear whether log index or index' 
correlates better with logao, but as the classical Mg-ao re- 
lation was established with Mg2 measured in mag we adopt 
this approach here for all other indices as well. The conver- 
sion between an index measured in A and magnitudes is 

index \ 



index' — —2.5 log ^1 



AA J (2) 

where AA is the width of the index bandpass (see e.g., W097 
and Trager et al. 1998 for a list of bandpass definitions). Fe3' 
is defined as 

Fe4383' + Fe5270' + Fe5335' 



Fe3' = 



(3) 



Fi gure [ic| shows index-ao relations for eight different 



log<Fe> (0.209 ±0.047) 
logH/3G (-0.081 ±0.042) 

logC24668 (0.429 ±0.096) 



(0.075 ± 0.025) 
(0.084 ± 0.042) 
(-0.231 ± 0.082) 
(-0.169 ± 0.038) 
(0.63 ±0.06) 



(1) 
(2) 
(1) 
(2) 
(1) 



Note — Errors are estimated by a Jack-Knife error analysis. 



metal i idices and two Balmer-line indices. The best fitting 
linear relations and the scatter are summarized in Table |^ for 
all indices considered in this paper. For the fits we used an 
ordinary l east square metho d, minimizing the residuals in y- 
direction (|lsobe et al. 1990|, hereafter OLS(Y|X)). Included 



in the fit are all galaxies with old stellar populations, i.e., 
all Es plus the bright SOs NGC1380 & NGC1381; in total 13 
galaxies. The 1-a scatter around the relation was robustly 
estimated by deriving a value which includes 9 out of 13 
galaxies (69%). A correlation coefficient derived from a (non 
parametric) Spearman rank-order test is given in the lower 
right corner of each panel in Figure 0. The probability that 
the parameters are not correlated is given in brackets. 

For the galaxies with old stellar populations the Mg- 
iTo relation is in excellent agreement with the literature 



References: (1) j0rgensen (1997, 11 nearby clusters) 
(2) j0rgensen (1999, Coma cluster) 

(J0rgensen 1997; CoUess et al. 1999). Remarkably, the Fe- 
line-indices also show a clear positive correlation with the 
central velocity dispersion and little scatter. This is the first 
time such strong correlations have been found at a signifi- 
cant level. We note that all the Fe-line and Ca4455-(Jo rela- 
■ tions show a slope consistent with a value of ~ 0.035. In con- 
trast the slope of the Mg-lines and C2 4668 are significantly 
steeper (see dot-dashed line in Figure ^ and Table ^ . 

Although our Mg2-(7o relation agrees well with the 
literature values we find significant differences for other 
log(mciea;)-(To relations compared to the data of j0rgensen 
(1997,1999). Table ^ shows a comparsion of the slopes. The 
log(<Fe>)-o-o relation seems to be far steeper in the Fornax 
cluster whereas the logH/^G-o-Q relation is shallower com- 
pared to Coma. The log C24668-(To relation in Fornax is 
marginally consistent with j0rgensen (1997). It is not clear 
why the \og{index)-ao relations for H/3g and <Fe> should 
be different to the Coma cluster. We will present a possi- 
ble explanantion at the end of this Section and in Section ^ 
where we discuss our results. 

The centres of the two bright and old SOs NGC1380 
& NGC1381 follow generally well the relation set by the 
elliptical galaxies. The lower luminosity SOs have velocity 
dispersions ao ^ 70 km s~^ and show a large scatter about 
the mean relation of the old galaxies. However, it is worth 
noting that they exhibit generally weak Mg absorption and 
some of the faint SOs show as much Fe absorption as L* 
ellipticals. Fornax A, the brightest galaxy in our sample has 
a central velocity dispersion of ctq ~ 220 km s~^ which is too 
low compared to ellipticals of this luminosity in the Faber- 
Jackson relation (see Figure ^l|) . It also departs significantly 
from the Mg-cro relations in the sense that it shows too weak 
Mg- absorption. As Fornax A is r egarded as the product of 
a recent merge r (Schweizer 198C; Schweizer 1981; Mackie & 



Fabbiano 1998) we interpret our results as strong indications 
of at least one young stellar component in this galaxy. 

One would expect the young stars in this galaxy to pro- 
duce strong Balmer absorption lines (as seen in Figure ^ 
and to dilute (or weaken) the metal lines of the underlying 
older stellar component. However, if the burst mass is not 
too small, the relative abundances of metal lines should in 
first order not be affected (see discussion of composite stel- 
lar populations in Section 5.3). Yet, we find that Fornax A 



deviates only from the Mg-ao relation and not from any of 
the other metal-index-cro relations (Figure p^). We interpret 
this as good evidence that the underlying older stellar popu- 
lation of Fornax A is significantly different from ellipticals at 
this velocity dispersion, i.e., the [Mg/Fe] ratio is lower, close 
to solar. 



© 2000 RAS, MNRAS 000, 



14 Harald Kuntschner 



0.40 ; 



0.30 r 



(a) 





0.10 



O.BO (0.10%) 




; 0.0 



X 



0.05 - 



0.00 [ 



-- -0.1 



-O.BO (0.10%) 



-0.2 



1.6 1.8 2.0 2.2 2.4 

log <^o /[km/s] 



1.6 1.8 2.0 2.2 2.4 
log /[km/s] 



2.6 



Figure 10. Selected metal absorption indices and Balmer line indices plotted against logfTQ. All indices are measured in magnitudes 
following the conversion from Equation ^ Note that all metal indices show a positive correlation and the Fe-indices show very similar 
slopes. The dot-dashed line indicates a slope of 0.035 centred at logcro = 2.2. Note that the Mg indices exhibit a much steeper slope. 
The Spearman rank-order correlation coefficient for each dataset is shown with the significance level in brackets in the lower right corner 
of each panel. All Es and the large SOs NGC1380 & NGC1381 were included in the fits. 
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Figure 11. The Faber-Jackson relation for the Fornax sample 
(assuming m — M = 31.2). The magnitudes are taken from the 
RC3 (de Vaucouleurs et al. 1991). The best Unear fit (OLS (Y|X) 
including all Es, NGC1380 & NGC1381) is shown as a solid Une. 



Two of the ellipticals stand out from the normal metal 
index-o-o relation: NGC1373 and IC2006 (labelled in Fig- 
ure |l^. These galaxies always show stronger metal line 
absorption than what would be expected from the mean 
relation. This is most prominent in the FeS'-ao diagram 
(panel c). We note however, that the galaxies follow the 
mean Faber-Jackson relation (Figure y_l|). There is little 
known about the galaxy NGC1373; perhaps the best expla- 
nation why this (elliptical) galaxy is somewhat off the mean 
relation is to regard it as a transition galaxy between the 
sequence of Es and the faint SOs. However, IC2006 has been 
studied in detail by Schweizer et al. (1989). They found a 
large counter-rotating ring of neutral hydrogen (HI) associ- 
ated with faint optical features and suggest that the HI ring 
may have formed during a merger which created IC2006. 
Franx et al. (1994) re-analysed the optical photometry of 
Schweizer et al. taking into account the inclination of the 
galaxy and concluded that it probably has a large disc in 
the outer parts which is seen almost face on and therefore 
difficult to detect. They suggest that it should be classified 
as E/SO rather than a bona fide elliptical. 

It seems plausible that the (perhaps peculiar) merger 
history of this galaxy is the reason for its deviation from the 
index-cTo relations. However from our data, it is not clear 
whether the stellar populations of IC2006 are too metal rich 
or whether the central velocity dispersion is reduced com- 
pared to other elliptical galaxies of this mass. If indeed this 
type of galaxy is more frequent in other clusters, such as the 
Coma cluster (see j0rgensen 1999), it would explain why 
previous authors did not find a clear correlation of Fe-lines 
with (70 . A detailed analysis of the kinematics and stellar 
population of this galaxy could be very valuable for our un- 
derstanding of how todays early- type galaxies were created. 

Panels (i) & (j) in Figure hd show the index-ao re- 



lations for two Balmer lines. Both indices show negative 
correlations. Elliptical galaxies and the bulges of NGC1380 
and NGC1381 show little spread around the mean relation 
whereas the younger galaxies, most remarkably NGC1316, 
tend to have significantly stronger Balmer absorption at a 
given (To. We emphasize here that the slope in the relation 
of the galaxies with old stellar populations is mainly caused 
by a metallicity effect (metal poorer galaxies have stronger 
Balmer absorption) and has little to do with age differences. 
The two "metal-rich" galaxies IC2006 and NGC1373 are de- 
viant from the main H7a'-o"o relation in the sense of lower 
H7A line strengths. This is caused by the residual metallicity 
sensitivity of H7A. The side-bands of this index are located 
on metal lines which lower the pseudo-continuum level and 
thus weaken the index. 



7 GLOBAL RELATIONS 

In this section we investigate the relations between our 
best age, metallicity, [Mg/Fel estimates and the central 
velocity dispersions. Figure |l^ presents the results. The 
ages and metallicities were estimated from a Fe3 vs H7A 
age/metallicity diagnostic diagram (Figure ^ in combina- 
tion with V96 models. The errors are evaluated following 
the procedure outlined in Section 5.1 but only including 
the Poisson error for individual galaxies. Some of the young 
galaxies are at the edge or outside the range of the model 
predictions which prevents an accurate error evaluation. For 
the latter galaxies we do not plot error bars. Notice that the 
ages and metallicities are derived parameters which carry all 
the caveats discussed in the previous sections. For example 
the independent measurement errors of the line-strength in- 
dices translate into correlated errors in the age - metallicity 
plane due to the residual age/metallicity degeneracy in the 
Fe3 vs H7A diagram. We note that the results presented in 
the following paragraphs would not change significanlty if 



H/3 is used as an age indicator (see Figure CI 



and 



C2 in 



the Appendix for a comparison of the age and metallicity 
estimates derived from H7A & H/3). The Mg-overabundance 
is estimated by evaluating the difference in metallicity esti- 
mate between a Mg2-H/3 and a Fe3-H/3 diagram (see Fig- 
ure K]3 in the Appendix for an estimation of [Mg/Fe] using 
H7aT 

Our age estimates of ellipticals do not show a signifi- 
cant correlation with log cto (panel a) . With the exception of 
Fornax A all galaxies with ao > 70 km s~^ show roughly the 
same age whereas the younger galaxies populate the low ve- 
locity dispersion range. However, there is a hint that the two 
dynamically hottest galaxies are younger than their smaller 
brethren. 

For galaxies with old stellar populations there is a clear 
correlation between the central metallicity and the central 
velocity dispersion ao (panel b) . Consistent with our findings 
for the Fe3-(Jo relation the galaxies IC2006 and NGC1373 
show a stronger metal content than what would be expected 
from the mean relation. The young SOs spread over the whole 
metallicity range. The best fitting OLS(Y|X) relation (solid 
line, Jack-Knife error analysis) to galaxies with old stellar 
populations gives: 



[Fc/H] = (0.56 ± 0.20) log (7 - (1.12 ± 0.46) . 



(4) 
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Figure 12. The global relations between log age, metallicity, [Mg/Fe] and logao are shown. The filled circles represent elliptical galaxies 
and the open circles stand for the SOs. The two open circles with a cross indicate NGC1380 and NGC1381. ES0358-G25 and ESO359-G02 
are represented by a triangle and star respectively. For galaxies which lie at the edge or outside the model predictions in the age/metall 
icity diagnostic diagrams we do not show error bars. A Spearman rank-order correlation coefficient is shown in the lower right corner of 
each panel (significance in brackets). The Spearman rank-order test includes all galaxies with old stellar populations. The dashed line in 
panel (b) shows a linear fit excluding NGC1373 and IC2006. The long-dashed line in panel (d) shows the relation from j0rgensen (1999) 
for the Coma cluster. Panel (g) shows the metallicity estimates derived from a H/3 vs Fe3 diagram (y-axis) plotted against the metallicity 
estimates derived from a H/3 vs Mg2 diagram (x-axis). 



A correlation coefficient derived from a Spearman rank- 
order test (including all ellipticals and the two large SOs) is 



given in the lower right corner of each panel in Figure 
The probability that the parameters are not correlatec 



12 



given in brackets. Excluding NGC1373 and IC2006 from the 
fit gives the following relation (dashed line in panel b): 

[Fe/H] = (0.82 ± 0.18) log a - (1.72 ± 0.40) . (5) 

In the age-metallicity plane (panel c) we find a statis- 
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tically significant relation in the sense that the more metal 
rich (and also more luminous) galaxies are younger. The 
slope of this relation is similar to what j0rgensen (1999) 
found for the Coma cluster (see also Worthey, Trager & 
Faber 1995), yet the Fornajc galajcies with velocity disper- 
sion (To > 70 km s~^ span a much smaller range in age and 
[Fe/H]. We note that the non-treatment of non-solar abun- 
dance ratios in combination with correlated errors could 
be the sole reason for the trend found in Fornax. The di- 
rection and magnitude of correlated errors for a galaxy of 
solar metallicity and 8 Gyrs age are shown in panel (c), 
top right corner. Following on from the age-met allicity re- 
lation j0rgensen (1999) established for the Coma cluster an 
age-[Mg/H]-(Jo relation. It would be very interesting to see 
whether such a correlation exists also in Fornax. However, 
the small number of galaxies combined with a rather small 
spread in age makes such an analysis very insecure and has 
therefore not been attempted. 



Th3 Mg-overabundance shows a weak positive corre- 
lation with central velocity dispersion and [Fe/H] in the 
sense that dynamically hotter and more metal rich galax- 
ies are more overabundant (panel d, f). In the Fornax clus- 
ter significant overabundances are found for galaxies with 
ao ^ 100 km s"^ or [Fe/H] ^ 0.0 (panel f, g). The best 
fitting linear relation between [Mg/Fe] and logCTo is: 



[Mg/Fe] = (0.49 ±0.18) logcr- (0.80 ± 0.41) . 



(6) 



This relation is qualitatively in agreement with the re- 
sults from the Coma cluster (j0rgensen 1999, long-dashed 
line in Figure p^). The scatter about the [Mg/Fe]-(To re- 
lation in Fornax is consistent with the errors for [Mg/Fe], 
but there seems to be a rather large spread in the [Mg/Fe] 
ratio at a given metallicity (panel f). Although the latter is 
in good agreement with our findings in the Fe-Mg2 diagram 
(Figure ^), we note that the errors are heavily correlated in 
the [Mg/Fe]-[Fe/H] diagram. There is no significant corre- 
lation of the [Mg/Fe] ratio with log age where the young SOs 
show solar or slightly less than solar [Mg/Fe] ratios. 



8 DISCUSSION 

In this study, great care was taken to calibrate the line- 
strength measurements to a standard system in which we 
can compare the results with theoretical model predictions 
(Section ^). The accuracy of this calibration is vital when 
one wants to derive absolute age and metallicity estimates. 
Although we were able to demonstrate the high quality of 
our calibration some unresolved issues such as the system- 
atic offset in the Mg2 vs Mg b diagram, the rather large 
rms error in the original Lick/IDS stellar library and per- 
haps most important of all the largely unknown effects of 
non-solar abundance ratios prevent us from deriving accu- 
rate absolute age and metallicity estimates. However, for the 
discussion of relative differences in the stellar populations of 
early-type galaxies our data set and current models are of 
excellent use. 

In this paper we have made use of two stellar popula- 
tion models provided by Worthey (1994) and Vazdekis et al. 
(1996). Both models make use of the Lick/IDS fitting func- 
tions but have otherwise somewhat different prescriptions to 



predict line-strength indices of integrated single-burst stel- 
lar populations (SSP). The predictions of the two models 
are consistent and our conclusions would not change if only 
one of them had been used for the analysis. To our knowl- 
edge, this would be also true if we had used any other model 
which makes use of the Lick/IDS fitting functions. 

One of the most important results from this study is the 
homogeneity of the stellar populations in dynamically hot 
early-type galaxies in the Fornax cluster. Apart from For- 
nax A all early- type galaxies (Es & SOs) with ao > 70 km s^^ 
are of roughly the same age and their central metallicity 
scales with logao- The homogeneity is refiected in tight re- 
lations of observables such as Mg-cro and Fe-ao and a clear 
correlation of [Fe/H] with the central velocity dispersion. 
The existence of the latter is reassuring in terms of our cur- 
rent understanding of the colour-magnitude-relation (CMR) 
in clusters being mainly a result of increasing metal l icity 



with increasing luminosity (Kodama & Arimoto 1997 



Fer- 



levi ch et al. 1999D 

Previous authors (iFishcr, Franx fc Illingworth 199^ : 



j0rgensen 1997; j0rgcnsen 1999) pointed out that the lack 
of a correlation of Fe-absorption strength with central veloc- 
ity dispersion would give evidence for a second parameter or 
conspiracy of age, metallicity and [Mg/Fe] ratio which keeps 
the CMR tight. For example in the Coma cluster j0rgensen 
(1999) did not find a strong correlation of <Fe> with central 
velocity dispersion and hence her [Fe/H]-(To relation is also 
not significant. However, both the <Mg>-(To and [Mg/H]-cro 
relation are clearly seen in Coma. In this context it is impor- 
tant to note that the slope of the Fe-ao relation which one 
would expect from the cha nge of metallicity in the CMR 
(Kodama & Arimoto 1997) is quite shallow and therefore 



only detectable with high S/N data. In contrast, the Mg-ao 
relation is steeper, and therefore easier to detect. The rea- 
son for this is a combination of a larger dynamical range in 
the Mg indices compared to the average Fe-index and an in- 
creasing Mg overabundance with central velocity dispersion 
giving a steeper slope than what would be expected from 
the change in metallicity only. 

An alternative explanation for the lack of a Fe-ao re- 
lation in Coma could be based on galaxies such as IC2006, 
which do not follow the Fe-o"o relation very well. If this type 
of galaxy is more frequent in the Coma cluster than in For- 
nax it would be impossible to find a clear Fe-ao relation. In 
summary we find that in the Fornax cluster there is no need 
for a second parameter such as age, metallicity or [Mg/Fe] 
to keep the CMR tight. Indeed, we favour an interpretation 
where small variations of age, metallicity and/or [Mg/Fe] at 
any given ao are responsible for some real scatter in the scal- 
ing relations for the Fornax cluster. However, we emphasize 
that this may not be true for other (larger?) clusters. 

In addition to the population of old, dynamically hot 
early-type galaxies, we find a sizeable fraction of young, dy- 
namically colder (ctq ^ 70 km s~^) systems within our mag- 
nitude limited survey. Some of the young SOs (NGC1375, 
ESO359-G02 and ES0358-G25) fit in remarkably well with 
the predictions of galaxy harassment in clusters (Moore, 
Lake, Katz 1998; Lake, Katz, Moore 1998). In this sce- 
nario, medium sized disc galaxies (Sc-type) fall into a cluster 
environment and get "harassed" by high speed encounters 
with cluster galaxies. The end-products are small spheroidal 
galaxies where some gas of the disc is driven into the centre 
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of the galaxy. This gas is likely to be turned into stars in 
a central stellar burst. We note, that most of these young 
galaxies are in the periphery of the Fornax cluster consistent 
with having been "accreted" onto the cluster from the field. 

Two of the SOs which show young populations in the 
centre, also have extended discs (NGC1380A and IC1963). 
This seems to be in contradiction with the harassment pic- 
ture. However we emphasize, that the existing harassment 
simulations do not include spirals with a substantial bulge 
component. Here the bulge is likely to stabilize the disc 
and the end-products may be able to keep substantial disc 
components (Ben Moore, private communication). The ex- 
istence of a population of dynamically colder galaxies with 
young stellar populations in the nuclear regions is in agree- 
ment with a typical (nearby) cluster CMR where one finds 
a tail of blue galaxies towards the faint end (e.g., Terlevich 
1998). Furthermore Terlevich et al. (1999) demonstrate for 
the Coma cluster, using line-strength analysis, that these 
blue galaxies contain young stellar populations rather than 
being metal poor. 

It seems that in the Fornax cluster significant amounts 
of young stellar populations are predominantly found in low 
luminosity (lenticular) systems. However, for a sample of 
Coma cluster early-type galaxies Mehlert (1998) found that 
relatively bright SOs spread over the whole range in age (Es, 
excluding the cDs, are found to be old). This of course raises 
the question whether morphology is the driving parame- 
ter for young stellar populations (only SOs are younger) or 
whether luminosity is the important parameter (low lumi- 
nosity E & SO galaxies are on average younger). Taking the 
results from Coma and Fornax together we would like to ar- 
gue that in clusters it is only the lenticular galaxies which 
show signs of recent star formation and that low luminosity 
lenticular systems are more likely to do so. The latter may be 
just caused by the recent accretion of these low luminosity 
systems onto the cluster. 

So far we have addressed the age and metallicity distri- 
butions in the Fornax cluster with the help of line-strength 
indices. However there is more detailed information on the 
star formation (SF) processes to be gained if one investigates 
the [Mg/Fe] abundance ratios. When new stars axe formed 
chemical enrichment is predominantly driven by the ejecta 
of SN la (main producer of Fe peak elements) and SN II 
(producing mainly alpha elements). However, SN la are de- 
layed compared to SN II which explode on short time-scales 
of < 10® — 10^ yr. Taking this into account there are mainly 
two mechanisms which determine the Mg/Fe ratio in galax- 
ies: (i) the star formation time scale and (ii) the fraction of 
high mass stars, i.e., the initial mass function (IMF) (sec e.g., 
Worthey et al. 1992) . As re-confirmed in this study, the ma- 
jority of cluster early-type galaxies show a trend of increas- 
ing [Mg/Fe] ratio with central velocity dispersion. Galaxies 
with young stellar populations and/or low luminosity galax- 
ies show roughly solar abundance ratios. Given that the most 
luminous galaxies are also the metal richest, we emphasize 
that any realistic star-formation models have to be able to 
produce metal rich and Mg-overabundant stars at the same 
time. 

In principle one can reproduce the observed trends of 
overabundances and metallicity with varying star-formation 
time-scales: large galax ies form within shorter time-scale s 
than smaller galaxies ( Bressan, Chiosi, fc Tantalo 1996 ). 



However, this leeds to extremely short star- formation time- 
scales for the most massive galaxies. A plausible way to solve 
this dilemma would be a varying IMF where massive galax- 
ies have a top heavy IMF and low luminosity galaxies show 
a more Salpeter like IMF. For further discussions of the mat- 
ter see also Peletier (1999) and Tantalo, Chiosi & Bressan 
(1998). 

Recently, Thomas & Kauffmann (1999, see also Thomas 
1999) presented preliminary results from their semi-analytic 
galaxy formation models for the distribution of [Mg/Fe] in 
galaxies as a function of luminosity. In this scenario lumi- 
nous ellipticals are the last to form and hence Thomas & 
Kauffmann find a trend that the [Mg/Fe] ratio decreases 
with increasing luminosity, opposite to the observed trend. 
In general it seems very difficult with the current stellar 
population models to reproduce the observed Magnesium 
strength, and therefore [Mg/Fe] values ([Mg/Fe] ~ 0.4) of 
luminous ellipticals (Greggio 1997, but also see Sansom & 
Proctor 1998). 



9 CONCLUSIONS 

We have measured the central line strength indices in a mag- 
nitude limited sample of early-type galaxies brighter than 
Mb = —17 in the Fornax cluster and have applied the mod- 
els of Worthey (1994), Worthey & Ottaviani (1997) and 
Vazdekis et al. (1996) to estimate their ages, metallicities 
and abundance ratios. We find that: 

(i) Elliptical Galaxies appear to be roughly coeval 
forming a sequence in metallicity varying roughly from 
—0.25 to 0.30 in [Fe/H]. This result is consistent with the 
conventional view of old, coeval elliptical galaxies where the 
metallicity scales with the luminosity of the galaxy. This is 
refiected in scaling relations such as Mg-ao. Remarkably, 
we could show that all other metal line-strength indices also 
clearly correlate with the central velocity dispersion. In fact 
all Fe-line-cro relations are consistent with having the same 
slope. 

(ii) Lenticular Galaxies have luminosity weighted 
metallicities spanning the whole range of SSP model predic- 
tions. Lower luminosity SOs show luminosity weighted ages 
younger than those of the ellipticals. However, the centres 
of the bright lenticular galaxies NGC1380 and NGC1381 re- 
semble the properties of ellipticals suggesting that they ex- 
perienced similar star formation histories. The peculiar SO 
galaxy Fornax A (NGC1316), which is the brightest galaxy 
in the sample, has strong Balmer lines implying a very young 
luminosity weighted age, yet the metallicity is equal to the 
most metal rich Es. This is consistent with Fornax A having 
been involved in a recent gaseous merger. The SOs NGC1380 
and NGC1381 follow the index-ao relations of the ellipticals 
very well. However, the SOs with a young stellar component 
generally show a large scatter around the scaling relations. 

(iii) Our conclusions are based on several age/metallicity 
diagnostic diagrams which give consistent results. Further- 
more we demonstrate the advantage of using an emission 
robust age indicator such as H7A when analysing the stellar 
populations of extragalactic objects. 

(iv) We have discovered that two of the fainter and very 
metal poor lenticular galaxies appear to have undergone ma- 
jor star-formation in the last 2 Gyrs (in one case very much 
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more recently). We note that, like Fornax A, most of the 
young galaxies lie on the periphery of the cluster. This is 
consistent with the harassment picture where these galax- 
ies are accreted from the field and undergo a morphological 
transformation with a central star burst. 

(v) The elhptical galaxies and the SO NGC1380 exhibit 
overabundances up to 0.4 dex in Magnesium compared to 
Fe. There is a trend that the most massive and metal rich 
galaxies are the most overabundant, whereas the fainter Es 
approach solar ratios. This trend is inconsistent with the 
currently available semi-analytical predictions for hierarchi- 
cal galaxy formation. SOs with young stellar populations are 
consistent with roughly solar abundance ratios, and may 
even be slightly underabundant. Remarkably also Fornax A, 
the brightest galaxy in our sample, shows close to solar abun- 
dance ratios which is not what one would expect of an early- 
type galaxy of its size. 

(vi) Furthermore we note that abundance ratio trends, 
which are not included in the models, can lead to a change 
of relative age and metallicity estimates depending on which 
index combination is used in the analysis. As long as the 
non-solar abundance ratios are not properly incorporated 
into the models the estimation of absolute ages of integrated 
stellar populations remains insecure. 
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APPENDIX A: KINEMATICS 



In Figure Al a literature comparison of the central velocity 



dispersion measurements is presented. The most recent data 
for Fornax galaxies from Graham et al. (1998) as well as the 
literature compilation of Smith (1998, Private communica- 
tion) is in good agreement with our data. The comparison 
with the literature compilation by McElroy (1995) shows a 
somewhat larger scatter. In order to establish an average 
error for the velocity dispersion estimates the mean scatter 
of the data compared to Smith (1998) and Graham et al. 
(1998) is evaluated. For the comparison all galaxies with 



(JO < 70 km s~^ and also NGC1419 (marked in Figure Al) 
from the Smith compilation were excluded. The mean scat- 
ter is 0.027 in logcro. Subtracting in quadrature a mean er- 
ror of 0.015 for the literature data gives an error of 0.022 [in 
logcro] for galaxies with ao > 70 km s~^. For galaxies with 
(Jo < 70 km s~^ we adopt the rms scatter of the template 
stars as estimate of the error (Alogao = 0.074). 
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Figure Al. Comparison of our central velocity dispersions with 
McElroy (1995, triangles). Smith (1998, squares) and Graham 
et al. (1998, diamonds). Error bars reflect the literature errors 
and the rms scatter due to different template stars for the AAT 
data. 



APPENDIX B: LICK/IDS CALIBRATION 



In Figure Bl a comparison between the original Lick/IDS 
index measurements of stars and selected galaxies in com- 
mon with our data is shown. The individual diagrams show 
the scatter (Lick/IDS - AAT) around the mean offset vs the 
average of Lick/IDS and AAT. 

Figure shows the velocity dispersion corrections for 
each index as derived from broadened stellar and selected 
galaxy spectra. 



APPENDIX C: COMPARISON OF DERIVED 
AGES, METALLICITIES AND [MG/FE] RATIOS 



Here we present three Figures (CI, C2 & C3) comparing our 
age, metallicity and [Mg/Fe] estimates derived with differ- 
ent index combinations. In general we find very consistent 
results when we use e.g., H/3 instead of H7A as an age indi- 
cator. Most of the differences, such as the [Mg/Fe] ratios of 
some of the young lenticular galaxies, is caused by the data 
points lying at the edge of the model predictions or where 
they have been extrapolated. The extrapolation was nec- 
essary because the models reflect roughly solar abundance 
ratios and cannot account for e.g., the strongest Mg2 ab- 
sorption found in our sample. 



APPENDIX D: FINAL CALIBRATED 
CENTRAL INDICES 

Dl Table of central index measurements 

The final, fully corrected central (2'.'x3'.'85) index measure- 
ments and associated errors for the Fornax galaxies and 
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Figure Bl. Comparison of 16 indices between Lick/IDS and our measurements after the Lick/IDS offset correction (see Table y) lias 
been applied. Open diamonds, filled circles and the filled triangle represent stars, Fornax-galaxies and NGC3379 respectively. The formal 
error in the offset is shown as an error bar in the upper right corner of each panel. 
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Figure B2. Velocity dispersion corrections. Dots and open circles represent stars and galaxies respectively. The solid line connects the 
mean of all data points in each c-bin. Note that for H/3 only stars with equivalent width of H/3 > 1.1 A have been used to determine the 
mean. See text for details. 
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Figure CI. Comparison of the ages derived from a Fe3 vs H7A 
diagram and a Fe3 vs H/3 diagram using the mod els of Vazdekis 
(1996). For details of the method see Section 5.1. 



Figure C3. Comparison of the [Mg/Fe] ratios derived from a 
Fe3 vs H7A, Mg2 vs H7A diagram and a Fe3 vs H/3, Mg2 vs H/3 
diagram using the rnodels of Vazdekis (1996). For details of the 
method see Section 5.1. 
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Table Dl. Wavelength definition of H/3g 
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Figure C2. Comparison of the metallicities derived from a Fe3 
vs H7A diagram and a Fe3 vs H/3 diagram using the mod els of 
Vazdekis (1996). For details of the method see Section 5.1. 



NGC3379 are presented in Table D2. For each galaxy we 
give our Lick/IDS index measurement in the first row and 
the lo" error in the second row. The second column (first 
part of Table D2) lists central velocity dispersions in log a 
units. The last column (part two of Table D2) lists H/3g (not 
Lick/IDS index); for further details see Section D2. 



Index 


Index bandpass 


Pseudocontinua 


H/3g 


4851.32 - 4871.32 


4815.00 4845.00 






4880.00 4930.00 



D2 H/3 vs H/?G relation 



Figure Dl shows a plot of H/3 equivalent width vs H/3g 
equivalent width. The bandpasses of the H/3g index are de- 
fined by j0rgensen (1997), based on an earlier definition of 
an H/3 emission index by Gonzalez (1993), in such a way 
that the influence of the Fe feature right next to the H/3 
absorptio n fe ature is minimized. For wavelength definitions 
see Table El . 



The indices H/3 & H/3g show an excellent corr elation for 
the Fornax sample. The solid line in Figure Dl shows the 
fit to all the data, excluding ES0358-G25 (open triangle) 
because of its emission contamination, using a method which 
bisects the ordinary least squares fits made by minimising 
the X and the Y residuals: 



H/3g = (0.862 ± 0.027) ■ H/3 + (0.568 ± 0.051) 



(Dl) 



The derived relation (errors from a Jack-Knife analysis) is in 
very good agreement with the relation found by J0rgensen 
(1997). 
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Table D2. Fully corrected Lick/IDS indices for the central 
(2'.'x3^'85) extraction 



Name log cr G4300 Fe4383 Ca4455 Fe4531 C4668 H/3 Fe5015 Mgi 
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Table D2 - continued 
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